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Abstract 

Let  T]  be  a  superprocess  under  the  measure  P.  We  show  the  existence 
of  probability  measures  which  are  absolutely  continuous  with  respect 
to  P,  and  whose  Radon-Nikodym  derivatives  are  suitably  normalized 
exponential  functions  of  the  self  intersection  local  time  of  77.  These 
measures  correspond  to  measure  valued  processes  exhibiting  a  cer¬ 
tain  amount  of  self  interaction.  A  finite  time  divergence  of  the  total 
mass  (1,7 ]t)  is  shown  to  occur  in  a  related  model  in  which  the  change 
of  measure  involves  the  occupation  measure  of  the  superprocess.  As 
an  independently  interesting  side  issue  we  also  obtain  a  number  of 
results  related  to  a  self-interacting  version  of  the  Feller  diffusion. 
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Introduction  and  main  results 


1.1  On  interactions 

Our  aim  in  this  paper  is  to  study  the  effect  of  certain  exponential  changes  of  measure 
for  inducing  an  internal  interaction  on  the  behaviour  of  superprocesses. 

The  basic  object  of  study  will  be  a  super  Brownian  motion  (SBM),  { t]t}o<t<T ,  with 
T  >  0,  and  finite.  SBM  arises  naturally  via  an  infinite  system  of  branching  particles 
that  alternately  perform  Brownian  motion  and  undergo  critical  branching,  for  a  de¬ 
tailed  description  of  which  we  refer  the  reader  to  Dawson  (1993)  and  the  references 
therein.  While  we  shall  be  motivated  by,  and  use  the  language  of,  this  particle  system 
(and  the  relation  between  additive  functionals  on  the  particles  and  the  superprocess 
as  in  Adler  1992)  we  shall,  throughout  the  paper,  work  exclusively  with  the  limiting 
superprocess. 

The  (self)  interaction  that  we  wish  to  apply  to  the  superprocess  is  one  which  will 
cause  the  “particles”  to  be  attracted  either  to  one  another,  or  to  a  particular  region 
of  space.  Broadly  speaking,  there  are  two  ways  in  which  to  do  this. 

The  first,  which  we  shall  not  develop  here,  but  is  interesting  in  order  to  understand 
some  of  our  results,  follows  from  the  fact  that  SBM  can  be  expressed  via  the  solution 
of  a  stochastic  partial  differential  equation,  describing  the  time  and  space  flow  of 
particles.  This  is  particularly  well  described  via  the  so-called  “historical  process”  of 
Dawson  and  Perkins  (1991)  and  Perkins  (1992). 

In  this  setting  it  is  possible  to  modulate  the  diffusions  of  individual  particles,  by 
having  them  attract,  or  repel,  one  another,  much  in  the  way  that  the  generator  of  a 
finite  system  of  independent  diffusions  is  adapted  to  yield  interaction.  This  approach 
was  developed  in  Perkins  (1992),  and  extended  to  more  singular  interactions  in  Adler 
and  Tribe  ( 1996a, b).  In  all  of  these  cases,  it  is  well  known  that  the  new  processes 
are  singular  (in  terms  of  the  probability  measures  they  induce  on  path  space)  with 
respect  to  the  original  SBM. 

The  second  approach,  which  we  adopt  in  the  current  paper,  involves  finding  more 
delicate  changes  that  leave  the  new  process  absolutely  continuous  with  respect  to 
the  old.  Attempts  to  do  this  via  extremely  delicate  changes  of  the  diffusion  seem 
to  be  doomed  to  failure,  (cf.  Adler  and  Ivanitskaya  1996  for  an  example.)  Thus  we 
adopt  an  approach  developed  initially  for  defining  a  self-avoiding  Brownian  motion 
(Varadhan  1969,  Westwater  1980,  1982).  The  idea  here  is  to  look  at  the  process  over 
some  fixed  time  interval  [0,T],  and  make  one  global  change  of  measure  that  favours 
interacting  paths.  Note  that  this  is  quite  different  from  the  above  approach,  since 
under  this  new  measure  the  new  process  will  no  longer  be  Markovian.  It  is  also 
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different  from  previous  changes  of  measure  dating  back  to  Dawson  (1978)  in  which 
only  the  branching  mechanism  of  the  superprocess  was  affected,  and  not  the  particle 
motions. 

To  describe  our  approach  more  precisely  we  need  some  notation. 


1.2  The  superprocess  and  the  Feller  diffusion 


With  T  >  0,  let  rj  be  a  canonical  process  on  the  space  D([0,T],  Mf($d)),  where 
denotes  the  space  of  finite  measures  on  5Rd.  For  any  m  £  and  a  real, 

let  Pm,T,a  be  a  probability  measure  on  D([0,T],  Mf($d))  such  that  under  Pm,T.a  the 
process  T]  is  Markov  and  its  distribution  is  characterised  by 


Em,T,ae-<L71t>  =  (1.1) 

where  Emji0l  denotes  the  expectation  under  Pm,T,on  f  £  B+,  the  space  of  non-negative, 
bounded,  measurable  functions  on  3£d,  and  u  is  the  unique  strong  solution  on  [0,  T] 
of  the  partial  differential  equation 

dut  _ 

dt 

u0  = 


(A  -f  a)ut  -  u\ ,  (1.2) 

/•  (1-3) 


Here  A  is  the  d-dimensional  Laplacian.  We  shall  write  P  and  E  for  PmiT,o  and  Emixto 
respectively. 


When  a  =  0  in  the  above,  the  resulting  measure  valued  process  is  super  Brownian 
motion  (SBM).  If  a  >  0  (a  <  0)  the  corresponding  process  is  SBM  with  additional, 
constant  rate,  mass  creation  (annihiliation).  The  reason  for  adding  this  extra  mass 
term  will  become  clear  later. 


We  denote  by  Xt  the  total  mass  of  the  process  r}t,  that  is 

Xt  -  (i ,m)  =  v($d)- 

As  is  well  known,  Xt  itself  is  a  nice  Markov  process  on  D([0,  T],  3£),  generally  known 
as  the  “Feller  diffusion”  and  satisfying  the  stochastic  differential  equation 

dXt  =  a  dt  +  y/xtdWu  (1.4) 

with  X0  =  (1,  rj0)  =  (1  ,m)  >  0,  and  W  a  standard  Brownian  motion. 

If  a  =  0,  then  X  dies  (i.e.  reaches  zero  and  stays  there)  in  finite  time.  If  a  >  0  then 
X  is  positive  for  all  t  >  0,  but  a.s.  finite  for  each  t.  If  a  <  0,  then,  since  can 
become  negative,  we  shall  stop  it  at  r0  :=  inf {7  :  Xt  =  0}. 

We  denote  by  PXtT,a  the  distribution  of  X  induced  by  Pm<T,a,  where  (l,m)  =  x  and 
denote  by  EXix,a  the  corresponding  expectation  operator. 
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1.3  Change  of  measure  via  SILT 

In  this  subsection,  set  a  =  0  in  (1.1),  so  that  we  are  dealing  with  a  standard  SBM. 
We  make  the  following  assumption  for  the  whole  of  this  subsection: 

Assumption:  In  addition  to  being  finite,  the  initial  measure  m  has  a  bounded  density 
with  respect  to  Lebesgue  measure. 

For  any  e  >  0,  an  approximate  self-intersection  local  time  (SILT)  for  the  superprocess 
■q  can  be  formally  defined  as 


Jr  =  fT  dt  [T  ds  l{|t-*|>£}  (d(x  -  y),  r]s(dx)rjt(dy )), 

Jo  J  0 

where  6  is  the  Dirac  delta  function.  J£  is  well  defined  for  d  <  7  (Dynkin,  1988). 
There  is  no  problem  sending  e  — »  0  in  dimensions  d  <  3,  and  renormalised  versions 
for  this  limit  exist  for  d  =  4,5.  (Adler  1992,  Rosen  1992). 

We  can  now  turn  to  our  first  change  of  measure.  Motivated  by  the  true  self-avoiding 
Brownian  motion  (or  “polymer  process”)  we  would  therefore  like  to  consider,  for 
d  <  7,  the  probability  measure  on  D({ 0,  T],  Mf(JRd)),  defined  by 


Qtm,T(dw )  — 


e^JT 


E{eXJr) 


P(dw),  A  >  0. 


Unfortunately,  however,  the  following  is  true: 


(1.5) 


Lemma  1.1  For  all  e,  T  >  0,  we  have  that  Em(eXJT )  =  oo. 


The  proof  of  the  Lemma,  which  depends  on  computing  the  moments  of  ,  is  defered, 
along  with  all  other  proofs,  to  the  following  section. 

As  a  consequence  of  Lemma  1.1,  there  is  no  way  to  work  directly  with  the  definition  in 
(1.5).  One  possible  way  around  this,  that  was  successfully  adopted  in  Adler  and  Iyer 
(1997)  to  construct  a  self-attracting  Brownian  motion,  is  to  consider  the  sequence  of 
measures 


<3"f  (*“)  =  j 


X”=o  Wt) 


E(  E”=0 


P(dw), 


n  >  1. 


(1.6) 


and  approach  the  measure  (1.5)  via  weak  convergence  as  n  — >  oo,  and,  perhaps,  also 
as  e  — >  0. 
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Dynkin  (1988)  has  shown  that  E(Jf)2  is  finite  for  d  <  7.  Conditions  for  the  existence 
of  the  second  moment  of  the  approximate  SILT  of  order  n  are  also  given  in  the  above 
paper.  Similar  moment  calculations  can  also  be  found  in  Rosen  (1992)  and  Feldman 
and  Iyer  (1996).  A  proof  of  the  fact  that  E^JfY  is  finite  for  d  <  7  and  all  j  >  1  is 
given  in  the  Appendix. 

The  main  result  of  this  paper,  which  indicates  that  a  program  as  described  above 
should  work,  is 

Theorem  1.2  The  collection  of  probability  measures  }e>o,n>i  is  tight  in  di¬ 

mensions  d  <  7. 


Thus,  given  {n^,  £k}k>u  the  sequence  of  measures  {Qm.T*’*}  contains  at  least  one 
subsequence  that  converges  weakly  in  £)([0,  T],  Mf(9ld)),  that  may  be  used  to  define 
a  limiting  polymer  measure.  We  shall  have  nothing  to  say  about  the  uniqueness  (or 
otherwise)  of  this  limit.  However,  we  shall  investigate  a  few  of  its  properties. 

Comment  for  Robert  We  may  have  something  to  say  about  uniqueness. 

Of  course,  while  Theorem  1.2  gives  us  a  candidate  for  a  new  polymer  process,  it  does 
not  tell  us  that  the  new  process  is  any  different  from  the  old.  (This  is  a  non-trivial 
issue,  in  view,  for  example,  of  the  aforementioned  problems  with  interesting  looking 
interactions  at  the  particle  level  that  disappear  in  the  infinite  density  limit.)  The 
following  result  indicates  that  the  polymer  measure  is  quite  different  from  the  initial 
superprocess  measure,  in  at  least  a  number  of  ways. 


Theorem  1.3  a.  For  all  A  >  0, 


iim  E'^(X,)  = 


40 


MO  +  ^  l  +  .uA'/'.m.lmll  )  ^  ^ 

|m|  4-  HzilzlZ  4  <  d  <7 


where  h\  and  hi  are  positive,  bounded  functions, 
b.  For  all  A  >  0  and  4  <  d  <  7 


lim  lim  E^xi^t)  =  Er, 


do 


(Xt) 


t  >  0. 


(1.7) 


(1.8) 


Note  firstly  that  for  the  basic  SBM  we  have  E(Xt)  =  |m|,  for  all  t  >  0,  and  all  d. 
Thus,  (1.7),  (1.8)  tell  us  that  not  only  is  the  polymer  measure  different  from  that  of 
SBM,  but  it  also  involves  “creation  of  mass”.  Looking  back  at  the  formulation  (1-6) 
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of  ,  and,  more  simply  if  not  so  precisely  (1.5),  it  is  clear  that,  under  Q,  SILT  is 
enhanced.  In  the  development  of  a  regular  polymer  measure  for  standard  Brownian 
motion  this  means  that  the  particle  path  must  be  changed.  Here,  however,  there  are 
two  ways  to  increase  SILT:  One  is  by  moving  the  paths,  the  other  is  by  affecting  the 
branching  so  that  the  total  mass  (and  so  SILT)  is  increased.  Theorem  1.3  indicates 
that  the  latter  is  certainly  happenning  here.  How  much  interference  with  the  particle 
motion  in  being  generated  by  the  change  of  measure,  and  how  the  two  effects  can  be 
uncoupled,  is  currently  unclear. 

Comment  for  Robert  Another  point  to  be  made  is  that  when  we  first  send  e  to 
zero  and  then  send  n  to  infinity,  the  expected  value  of  the  total  mass  of  the  limiting 
process  has  a  jump  at  t  =  0,  and  then  it  stays  constant  as  opposed  to  the  case  of  fixed 
n.,  when  there  is  a  continuous  increase  in  the  expected  total  mass.  This  suggests  that 
things  may  be  different  in  the  limit.  This  is  similar/a  contrast  to  the  fact  that  even  if 
we  start  with  an  absolute  continuous  mass  the  SBM  immediately  becomes  singular. 

A  second  rather  interesting  consequence  of  (1.7)  is  that,  at  least  for  dimensions  4  < 
d  <  7,  we  do  not  recover  the  basic  SBM  by  sending  the  interaction  strength  parameter, 
A,  to  0. 

As  an  interesting  aside,  note  that  even  the  case  of  fixed  n  and  e  =  0  gives  a  measure 
different  from  that  of  SBM.  This  result  is  at  variance  with  the  result  for  the  case  of 
the  attractive  polymer  model  for  the  Brownian  motion  (admittedly  a  different  -  but 
similar  looking  -  problem)  where  keeping  n  fixed  and  sending  e  to  zero  lead  us  back 
to  the  Brownian  motion  (Adler  and  Iyer  1997). 

Theorem  1.4  The  underlying  Feller  (total  mass)  process  is  uniquely  determined  un¬ 
der  the  limiting  probability  measure  lim^oo  limc|0  Qm^Ti  for  T  >  0  and  A  >  0  in 
dimensions  4  <  d  <  7. 

Despite  the  fact  that  Jf  is  not  exponentially  integrable,  we  do  not  lose  the  exponential 
integrability  of  the  basic  superprocess  in  the  following  particular  case. 

Theorem  1.5  Let  4  <  d  <  7  and  f  £  B+.  Then  the  limits  of  the  sequence  E^t  (e</ ,m'>) 
are  finite  when  we  first  let  e  0  and  then  take  the  limit  n  — >  oo,  for  t  <  (C'||/||0o)~\ 
where  C  =  1  +  2(1  V  H"1)^  V  T). 

1.4  Change  of  measure  via  occupation  measure 

Given  the  fact  that  when  changing  the  probability  via  SILT  we  (inadvertently?)  in¬ 
troduced  a  mass  creation  term,  it  seems  interesting  to  look  at  a  simpler  model  as  well. 
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Thus,  in  this  section  we  shall  introduce  interactions  in  such  a  way  that  the  process 
favours  more  mass  in  certain  locations  of  the  space. 

For  <j>  €  B+ ,  define  the  occupation  measure  (or  “man  hours”  cf.  Iscoe  1986a, b)  process 

Ft(<f>)  =  :=  [  {<f>,Vs)ds  (1.9) 

Jo 

and  consider  the  new  measure  defined  by 

~  \  jt 

QlUM  ■■=  A  >  0.  (U0) 

If,  for  example,  <f)  has  compact  support  A,  then  what  we  have  defined  is  a  measure 
valued  process  that  favours  realisations  putting  most  of  their  mass  in  A.  Again,  what 
is  not  clear  is  whether  this  is  done  by  redirecting  the  individual  particles  to  drift 
towards  A ,  or  by  generating  mass  creation  in  A.  In  SILT  model  we  saw  an  increase 
in  total  mass,  but  could  not  say  much  beyond  that. 

However,  if  we  take  <j>  =  1  then  some  simple  analysis  of  the  total  mass  process  (i.e. 
the  Feller  diffusion)  shows  that  the  mass  creation  is  so  intense  that  in  general  the 
process  will  explode  in  fixed  time.  This  is  described  below  in  the  following  section, 
which  is  essentially  independent  of  the  remainder  of  the  paper. 

A  more  careful  analysis  of  general  <f>,  and,  in  particular,  of  the  particularly  interesting 
case  when  <f>  =  5,  so  that  the  occupation  time  is  replaced  by  the  local  time,  will  be 
the  object  of  a  future  study. 


1.5  The  reinforced  Feller  diffusion 

As  before,  let  Xt  =  (l,Vt)  be  the  total  mass  process:  i.e.  the  Feller  diffusion  (1.3)  and 
fit  :=  Jq  Xsds  —  fit(  1).  Consider  the  following  change  of  measure,  which,  along  with 
the  remainder  of  this  subsection,  has  nothing  to  do  with  superprocesses: 

X  >  °-  (U1) 

Since  this  change  of  measure  clearly  defines  a  process  that  likes  to  take  high  values 
we  shall  call  it  a  reinforced  Feller  diffusion.  (Of  course,  the  term  “diffusion”  here  is 
somewhat  inappropriate,  since  this  process  is  not  Markovian!) 

It  is  actually  possible  to  explicitly  compute  the  Laplace  transform  of  Xt  for  0  <  t  <  T, 
under  the  new  measure,  and  the  result  is  surprisingly  dependent  on  the  values  of  a 
and  A.  Here  are  some  sample  results: 
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Theorem  1.6  Let  A,c  >  0,  and  set  7  =  A  —  a2  /  4.  Then ,  for  all  0  <  t  <T } 

_  exp{VA  stan^/Xt  -  -‘^)|} 

£*'r'o(e  1  -  expfv^xtanfVAT)}  1  ’ 

(ii)  If  7  >  0,  then 

f  r-  ,  r,  _1  /c-'/ytan(tan-1(5^.)+\/'r(T-0)  ^ 

exp{— ,/y  x  tanftan  *( - Z; - yV'U 

Fx  (e-cXt\  _  F1  _ _ _ _ _ ^5 - L—iL  (1.13) 

-&x  X  <y\^  )  r  i —  .  r.  —  1  /  rv  \  .  / — ml  1  V i  J } 


exP{\/T  a:  tan  [tan  +  \/^T}} 


(in)  If  7  =  0,  then 


,  c(l-f(T-t))-f  i 

exPl  l-f(T-t)+<%(l-f(T-i))-f) 


EXI,TA‘-eX‘)  = 
(iv)  If  y  <  0,  let  (3  =  Then 


(T— t)+<(c(l— f(T— «))- 
a 

exp{i_2or  a;} 


(1.14) 


r  c+w?._t  ■, 

pa  ,  -cxn  =  exp{— •W;  a:} 

ar,T’0^  exp{— X}  ' 


where, 


<-!  =  * 


'/9(e»"-l)  +  (c-f)(ew‘+l) 
[/3(eM'  +  l)  +  (c-f)(e2's,-l) 


(U5) 


(1.16) 


An  unexpected  corollary  of  the  above  is  that,  under  the  new  measure,  if  T  approaches 
a  particular  (a,  A  dependent)  time,  mass  creation  is  so  fierce  that  the  corresponding 
processes  diverge  to  +oo  in  finite  time.  More  precisely, 

Corollary  1.7  (Explosion  of  mass)  Let  7  be  as  in  the  previous  theorem.  Let  A  >  0. 

(i)  Ex,t, ,„(e-*)  0  as  T  |  ;Jj. 

(ii)  If  7  >  0,  then  ElTJe-‘x-)  -*  0  as  T  T  (f  -  tan'1^))^. 

(mj  If  a  >  0  and  7  =  0,  then  E^Ta(e~cXt)  — ►  0  as  T  t  2/ar. 

(t»j  If  7  <  0  and  a  >  2/3,  then  EfTa{e~cXt)  — *  0  as  T  f  ^ln  • 

A  number  of  additional  properties  of  the  reinforced  Feller  diffusion  are  given  in  Sec¬ 
tion  3  below.  While  they  seem  to  be  of  independent  interest,  and  so  are  included 
there,  they  are  not  related  to  the  main  theme  of  this  paper. 
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2  Proofs  and  supporting  results 


Given  the  form  of  our  change  of  measure,  it  should  come  as  no  surprise  that  most  of 
the  proofs  of  the  results  in  subsection  1.3  will  be  based  on  an  analysis  of  the  moments 
of  SBM. 

Dynkin  (1988)  derived  what  is  now  a  well  known  explicit  representation  for  these 
moments  in  terms  of  the  transition  densities  pt(x,y)  and  binary  graphs.  Since  his 
representation  will  form  the  backbone  of  the  proofs,  we  include  it  here  for  the  sake  of 
completeness. 

2.1  Moment  Formulae  for  the  Superprocess 

Let  T>k  be  the  set  of  all  directed  binary  graphs  with  k  exits  marked  1, 2, . . . ,  k.  Given 
such  a  graph  G,  let  L  be  the  set  of  directed  links  and  V  be  the  set  of  all  vertices.  If 
the  link  /  £  L  goes  from  vertex  v  €  V  to  vertex  w  G  V,  we  write  v  =  i{l),w  =  f(l). 
We  associate  with  each  element  v  6  V  two  variables, 

{sv>yv )  £  x 

which  we  refer  to  as  the  time  and  space  coordinates  of  v.  We  denote  by  V-  the  set  of 
all  entrances  of  G;  i.e.  the  set  of  all  vertices  v  (E  V  such  that  no  arrow  ends  at  v  and 
there  is  only  one  arrow  that  begins  at  v.  For  any  v  €  V-  we  set  sv  =  0  and  yv  =  xv. 
If  v  is  the  exit  labeled  by  j,  1  <  j  <  k  —  that  is,  v  is  such  that  there  is  exactly  one 
arrow  which  ends  at  v  and  no  arrows  begin  at  v  -  we  set 

(•s  vi  yv)  = 

Finally,  let  V0  denote  the  set  of  internal  vertices;  i.e.  those  vertices  that  are  neither 
entrances  nor  exits. 

Let  pt(x,y)  =  Pt(x  —  y)  =  (27ri)-d/2exp(— ||x  —  j/||2/2),  t  >  0,  be  the  standard 
Brownian  transition  density,  and  set  pt  =  0  if  t  <  0. 

Theorem  2.1  Let  p  be  a  super  Brownian  motion.  Then,  for  all  k  >  1,  t\, . . . ,  >  0 

and  €  B+, 

£(n  (/..’)«.>)  =  e  oD,  (2i) 

\i= 1  /  Devk 

where 

f  k 

Cd  =  /  n  m(dxv)I[Psm-si(i)(yf(l)  ~yw)  n  dsvdyvY[fi(zi)dzi.  (2.2) 

J  v€V_  leL  veVo  i=  1 
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For  detailed  examples  of  how  these  formulae  are  applied  see,  for  example,  the  original 
paper  of  Dynkin  (1988)  and  Rosen  (1992).  Adler  (1992)  explains  the  relationship  of 
Theorem  2.1  to  the  particle  picture. 


2.2  Preliminary  (but  crucial)  estimates 

In  most  of  the  calculations  that  follow,  we  shall  be  interested  in  deriving  moment 
estimates  in  which  the  crucial  issue  will  be  identifying  (negative)  powers  of  the  pa¬ 
rameter  t.  It  can  be  seen  from  the  proof  of  Theorem  1.3  that  the  exact  computation 
of  the  moments  are  both  extremely  involved,  and,  fortunately,  unnecessary  for  our 
purposes.  Since  most  of  the  techniques  involved  in  estimating  these  moments  have 
already  been  explored  in  great  detail  in  Dynkin  (1988)  and  Rosen  (1992),  the  proofs 
tend  to  be  a  bit  sketchy  at  times,  although  we  do  hope  that  we  have  covered  all  the 
important  points. 

We  start  with  some  general  moment  bounds,  in  which  /(e)  ~  g(e)  means  that  there 
exists  a  constant  c  such  that  lime_>o  /(e)/^(e)  =  c- 

Lemma  2.2  For  d  <  7, 

EMtY  ~  E  C(m,d,D)(lnT)^r-^+a(D’d)6(d,P,e),  (2.3) 


where  0  <  k(D)  <  j ,  0  <  c(D)  <  j —  k(D),  and 


a(D,d)  =  | 

r  \d\ 

{  \D\  -  2 HD) 

1  <  d  <  3 

4  <  d  <  7 

(2.4) 

f  1 

1  <  d  <  3 

b(d,D,e)=< 

(ln(e))‘<D> 

d  =  4 

(2.5) 

1  e-(d-4)k(D)/2 

5  <  d  <  7 

'y(D)  is  non-negative  if  d  is  even  and  zero  otherwise.  While  we  will  indicate  were  this 
extra  logarithmic  term  in  (2.3)  arises  when  d  is  even  in  the  proof,  we  shall  avoid  the 
extra  algebra  involved  by  treating,  in  detail,  only  the  cases  when  d  is  odd. 

Proof  of  Lemma  2.2.  Recall  that 

(j£)j  =  jQ  /  (%  -  y\ (2-6) 

i— 1  J 
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(2.7) 


Em(jy=  e  c'd, 

Dev7 } 

where, 

•?'  rT  /-T  y  7 

cd=h  /  iit2i-t2i-ii>e^2l^2i-i  /  n  -(^)n^(o-,<„(^(o-^o)  n 

i=r°  JO  j  veV_  leL  «€v0  i= 1 

(2.8) 

where  on  account  of  the  delta  function  we  have  yj m  =  Z{  if  sjm  =  t2i,t2i-\. 

Note  first  that  for  the  graphs  that  contribute  to  E(Jcy  the  exit  nodes,  in  the  notation 
of  (2.2)  are  all  labelled  (tk, z;),  k  —  1, . . .  ,2 j,  i  =  1, . . . ,  j. 

We  say  that  two  exit  nodes  ( )  and  (te,  zj)  are  paired  if  there  is  an  internal  node 
labelled  ( s,y )  such  that  there  is  an  edge  connecting  ( s,y )  to  (4,2;)  and  (t/,Zj).  Here 
is  a  preliminary  calculation  for  paired  nodes: 

Consider  the  term  corresponding  to  any  graph  D.  Let  the  graph  D  have  1  <  r  <  2 j 
entrance  nodes;  i.e.  r  sub-trees.  Let  |D|  denote  the  number  of  edges  in  D.  Suppose 
for  some  k  =  1  ,...,j  the  nodes  (t2k- \,Zk)  an<i  (^2 fc,^)  are  paired.  (Note  that  the 
space  variables  associated  with  the  exit  nodes  with  time  variables  t2k- 1  and  t2k  are 
identical  as  a  result  of  the  delta  function  appearing  in  the  definition  of  Je.)  We  shall 
call  nodes  paired  this  way  as  twin  pairs.  For  twin  pairs,  integrating  with  respect  to 
Zk  (and  applying  the  Markovian  nature  of  pt)  leaves  the  external  integral  involving 
f 2k— 1  and  t2k  in  the  form 

(27r)-<*/2'  L  L  +t2k-2s)~d/2  =c1(d,T-s,e)  (2.9) 

If  d  =  4,  then  ci(d,T  —  s,  e)  ~  ln(e)/2.  If  5  <  d  <  7,  then  cj(d,  T  —  s,  e)  ~ 
(4 /(d  —  4)(d  —  2))e-C-4^2.  Thus,  if  4  <  d  <  7,  (2.9)  diverges  as  e  —*  0  and  the 
important  terms  are  those  that  diverge  in  e.  If  1  <  d  <  3,  then  (2.9)  is  bounded  in  e 
and  we  can  take  e  to  be  zero,  to  get  cj(d,  T  —  s,e)  «  (2^_6^2/(d  — 4))(T  —  s)~C-4)/2. 

With  this  preparation  behind  us,  we  start  the  main  computation  for  evaluating  a 
typical  term  in  (2.7): 

Integrate  out  all  the  j—  space  variables  corresponding  to  the  2 j  exit  nodes.  For  all 
exit  nodes  that  are  twin  pairs,  also  integrate  with  respect  to  the  time  variables  as  in 
(2.9).  Suppose  there  are  0  <  k(D )  <  j  twin  pairs. 

Integrate  out  all  the  internal  space  variables  as  shown  in  the  appendix.  This  is 
standard  fare.  Ultimately  we  arrive  back  to  the  outermost  integral  with  respect  to 
m{dx ),  and  here  we  have  to  work  a  little  harder,  since  while  in  the  appendix  (as  is 
the  usual  case  for  these  computations)  we  are  interested  in  general  upper  bounds  here 
we  shall  require  more  precision. 
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Recall  the  assumption  that  m(dx)  =  f(x)dx,  for  some  /  €  B+(l Rd).  A  generic  outer¬ 
most  term  will  now  look  like 

J  m(dx)ps(x,y)gD(s,y)dsdy  =  j  Tsf(y)gD(s,y)dsdy ,  (2.10) 

where  Tt  is  the  heat  semigroup  and  go{s,y)  is  the  term  obtained  by  integrating  out 
all  the  spa.ce  and  time  variables,  except  for  the  space/time  variables  corresponding  to 
the  first  edge  from  the  root  in  the  sub-tree  with  entrance  node  labelled  (0,  a:). 

^From  the  appendix  we  know  that  /  gD(s,  y)dsdy  <  oo.  Moreover  gD(s,  y)  and  Tsf(y) 
are  strictly  non-negative  functions  on  [0,T]  x  $d  and  (0 ,T]  x  respectively.  gD  and 
Tsf  are  continuous  functions  on  (0,  T}x$d.  Let  I\n  be  any  sequence  of  compact  sets 
increasing  to  Then,  for  some  sequence  (sn,j/n)  (2.10)  is  equal  to 

.  lim  / Tsf(y)gD{s,y)lKn(y)dsdy  =  lim  TSnf(yn )  f  gD(s,y)dsdy ,  (2.11) 

n— >oo  J  n— +oo  j 

a  fact  which  follows  from  the  mean  value  theorem  applied  to  the  function  Tsf(y ) 
integrated  against  the  measure  gD{s,y)li<n(y)dsdy.  Since  0  <  TSnf(yn)  <  ||/||oo  we 
have  that  TSnf(yn )  converges  to  some  cp  €  (0,  ||/||oo)  as  n  — ^  oo  (since  the  right  side 
of  (2.11)  converges  to  the  left  side  of  (2.10)).  Consequently, 

m(dx)ps(x,y)gD(s,y)dsdy  =  cD  /  gD(s,y)dsdy. 

J  a[0,T]x5ftci 

If  we  use  the  above  technique  judiciously,  as  in  the  appendix,  we  can  integrate  out 
all  the  transition  density  functions  associated  with  the  internal  edges  while  integrat¬ 
ing  out  the  internal  space  variables.  Now  we  can  integrate  out  the  remaining  time 
variables  corresponding  to  the  remaining  exit  nodes  in  a  manner  similar  to  (2.9).  In 
evaluating  these  integrals  we  can  take  e  =  0.  If  s  =  (s/t, . . .  ,Sq._fc(D))  are  the  time 
variables  corresponding  to  the  internal  vertices  of  the  graph  D  that  are  connected  to 
these  exit  nodes,  we  can  write  this  factor  as 

n 

c3(D,d) 

where  r,-’ s  are  linear  combinations  of  T  and  the  components  of  s,  and  X)  ei  —  —c(D)d/ 2+ 
2(j  —  fc(jD)),  and  0  <  c(D)  <  (j  —  k(D))/ 2.  Note  that  each  will  occur  as  an  argu¬ 
ment  in  exactly  one  function  r,-.  In  case  d  is  even,  it  is  here  that  we  will  get  factors 
containing  logarithms  of  rt-’s.  This  will  make  the  subsequent  calculations  very  messy. 

A  note  here  on  the  exact  value  of  c(D).  This  will  be  used  later  in  the  proof  of 
Theorem  1.2.  To  compute  c(D),  we  integrate  as  described  above  all  the  twin  pairs. 
We  can  also  integrate  out  the  space  variables  corresponding  to  the  exit  nodes  for 
which  one  of  them  is  on  a  subgraph  with  only  one  edge.  Now  consider  the  residual 
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graph.  In  this  graph,  suppose  that  the  exit  nodes  are  labeled  z\,. . .  ,zp.  We  partition 
z}, ...  ,zp  into  disjoint  sets  5j , . . . ,  S',  such  that  the  exit  nodes  in  the  distinct  sets  a.re 
on  different  subgraphs.  Then,  c(D)  =  q.  Since  of  the  remaining  j  —  k(D)  pairs  of 
nodes  there  are  no  twin  pairs,  the  maximum  value  of  C(D)  will  be  as  given  above.  For 
example,  for  the  graph  in  figure  5  in  appendix  c(D)  —  1.  For  the  graphs  (i),(ii)  and 
(w)  in  figure  4,  c(D )  =  0,  whereas  for  the  graphs  (Hi),  (in),  (vi)  and  ( vii )  c(D)  =  1. 

This  leaves  us  with  \D\  —  2 j  internal  time  variables  that  on  integration  yield  a  factor 
T\°\~2j /c4(D).  For  example,  in  the  case  r  =  2 j,  04(D)  =  1.  If  r  =  j  and  all  the 
exit  nodes  are  twin  pairs,  then  there  are  j  internal  time  variables  and  04(D)  =  c? . 
Note  that  the  exponent  of  T  in  our  estimate  will  be  —c(D)d/ 2  +  \D\  —  2 k(D),  for 
4  <  d  <  7,  and  will  equal  —c(D)d/ 2  +  \D\  for  1  <  d  <  3.  In  case  d  is  even,  there  will 
be  additional  terms  in  ln(T). 

We  now  claim,  and  will  prove  below,  that 

j  <\D\-2k(D)<4j.  (2.12) 

It  is  now  easy  to  check  that  the  bounds  of  the  Lemma  are  obtained  by  simply  putting 
e  =  0  in  all  the  terms  that  do  not  diverge  as  e  — ►  0.  In  fact,  we  can  now  write  the 
following  estimate  for  Em(Jey  : 

Em(JJ  «  ]T  G(m,d,D)T-^+a^b(d,D,e),  (2.13) 

Dev  2j 


where  b(d,D,e)  equals  1  if  1  <  d  <  3,  equals  (ln(e))fc^^  if  d  =  4  and  is  given  by 
e-(d-4)k(D)/2  if  5  <  d  <  7.  k(D)  is  the  number  of  twin  pairs  in  D.  a(D,  d)  is  given  by 
\D\  if  d  =  1,3  and  equals  \D\  —  2 k(D)  if  4  <  d  <  7,  and  0  <  c(D)  <  j. 

This  proves  Lemma  2.2,  except  for  the  inequality  (2.12). 

Since  \D\  —  2k  <  |D|  <4 j,  the  right  hand  inequality  of  (2.12)  is  trivial.  For  the 
left  hand  inequality,  observe  that  we  can  associate  atleast  three  distinct  edges  of  D 
with  each  twin  pair,  and  that  there  are  k  twin  pairs.  Furthermore,  with  each  of  the 
remaining  exit  nodes  numbering  2 j  —  2k  we  can  associate  atleast  one  distinct  edge. 

Thus,  subtracting  2k  from  |Z>|  corresponds  to  taking  into  account  only  one  of  the 
three  edges  associated  with  the  paired  nodes,  so  that 

\D\  -2k>  (3  k  -  2k)  +  (2 j  -  2k)  =  2 j  -k>  j 

This  proves  the  claim.  ® 
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2.3  Proofs  of  the  main  results 

Proof  of  Lemma  1.1.  Since  J £  is  positive,  it  is  enough  to  bound  from  below  the 
contribution  to  .  We  do  this  by  restricting  our  attention  to  complete  binary 

graphs.  The  case  j  =  2  of  such  a  complete  binary  graph  is  illustrated  in  (2.3).  There 


Figure  1:  Complete  binary  graph  with  4  exits  (j=2) 


will  be  2 j\  terms  associated  with  a  complete  binary  graph.  This,  together  with  the 
estimate  for  obtained  above,  gives  us  that 

%Em,T(Jiy  >  ^C(m,d,D)b(i,D,e)(T-i'wy(2j)\  =  0(f),  (2.14) 

J'-  J- 

for  large  j.  This  implies  that 

-*  oo,  as  n-roo.  (2.15) 

j= o  *7* 

This  proves  Lemma  LI.  ® 


Proof  of  Theorem  1.2.  In  order  to  prove  tightness  of  the  family  of  measures 
WmlrhnM.oo,  we  start  with  a  bound  for  E^rdf ,Vt)k)-  Recall  that 


ZU  fEjiLm ^yyyy) 
E“=i  Em(J*Y 


(2.16) 


We  will  compare  Em((f,riT}k(J^)j )  with  £;m(«/f)J;  be.  the  corresponding  terms  in  the 
numerator  and  the  denominator.  Clearly, 


£»((/,  <  ||/||‘S.«l.^)‘W)i)  (2.17) 
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Let 


(2.18) 


Em(JT)j  -  53  Cd\  ar>d  Em((l,T]t)k(J^y)  —  ^3  Cd'i 

Dev2j  D'ev2]+k 

where,  as  before,  V2 j  is  the  set  of  all  binary  graphs  with  2 j  labelled  exits. 

Any  D'  €  V2 j+k  can  be  obtained  from  a  unique  D  G  V2j  in  one  of  the  following  two 
ways: 

(1)  Add  a  graph  D"  G  Vk  to  the  graph  D  as  a  seperate  component.  The  sum  over  all 
graphs  D'  G  V2j+k  obtained  in  this  way  will  give  a  term  of  the  form 

Em(l,m)kEm(J^y.  (2.19) 

(2)  Take  a  graph  D"  G  Vk .  Let  D"  have  1  <  r  <  k  components  (that  is  r  entrance 
nodes).  We  attach  1  <  p  <  r  of  the  components  of  D"  to  the  edges  of  D.  This  can 
be  done  in 


pi\  a(a  -  !)•  •  •  (a  -  1+  1) 


(2.20) 

ways,  where,  A (p,l)  =  {pi,...,pi  :  pi  >  0,i  =  l,..-,/,E|=i Pi  =  p}  and  a  =  \D\. 
Truly  speaking,  (2.20)  is  an  over  estimate,  since  one  has  to  divide  in  the  innermost 
sum  by  the  factorials  of  the  number  of  pi  s  that  are  identical. 


Observe  that  the  space  variables  in  the  inserted  graph  D"  can  be  integrated  out,  since 
the  function  associated  with  the  exit  nodes  in  this  graph  is  the  constant  function  1. 
After  doing  this,  the  space  variables  of  the  p  internal  nodes  at  the  points  where  the  p 
components  of  D"  have  been  appended  to  the  edges  of  D  can  also  be  integrated  out. 
Next  we  deal  with  the  time  integrals  coming  from  the  inserted  part.  If  D"  and  D  are 
as  illustrated  in  Figure  2.  then,  inserting  D"  into  D  on  the  edge  going  from  the  node 


(0,  x,) 


D 


// 


(  0,  x2) 


Figure  2:  Two  graphs  pre-insertion 
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(0,  x,  ) 


(t,  Zj)  (t,  z 2) 


Figure  3:  Graphs  of  Figure  2,  post-insertion 

labelled  («i,  2/1)  to  the  node  (52,2/2)  gives  the  graph  illustrated  in  Figure  3.  Note  that 
the  node  (0,®)  has  been  changed  to  (s4,x  1).  The  term  corresponding  to  the  graph 
D"  (considered  as  a  seperate  graph)  is 

J  m(dx)  J  ds  j  dy  ps(x,  y)  J  pt_s(y,zi)pt_s(y,z2)dzi<iz2  =\  m  \  ds.  (2.21) 

Thus  the  factor  corresponding  to  the  inserted  part  D"  into  D  as  shown  in  Figure  2.3 

will  be  t  t 

I  ds<  f  ds.  (2.22) 

Js4  Jo 

Consequently,  replacing  the  lower  bounds  in  all  the  time  integrals  of  the  first  internal 
node  from  the  root  of  any  component  of  D"  by  zero  (in  the  above  example  we  replace 
the  lower  bound  of  the  integral  w.r.t.  s  from  S4  to  0)  we  get  a  factor  that  is  bounded 
by  \m\-pEm(l,rit)k  \D",  where  Em(l,pt)k  \D»  is  the  contribution  to  Em(l,rjt)k  coming 
from  the  graph  D" .  This  leaves  us  with  a  graph  that  looks  identical  to  D  except  for 
an  additional  p  time  integrals  (created  at  the  points  of  insertion  of  the  p  components 
of  D"  into  D ). 

Carrying  out  the  estimation  of  the  term  as  in  the  proof  of  Lemma  2.2  and  using  (2.20) 
we  obtain 

Em(i,pt)k(jy  <Em(i,Vt)kEm(jy+  £  E 

Dev2j  D"evk  P=i  \  r  / 
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Pi  *  *  •  -  P/I  X 


*££ 

/=1  A  (p,/) 


P  \  P~Pl 


Pi 


P2 


p  -  ...  -  Pl—2 

Pi- 2 


'T'h+P 


xg  c(W.fl,«.0?ri-1)(4  +  2)...(k  +  rt 


,  (2.23) 


m  r-p ' 


where  I  =  I(D ,  jD",pi, . . .  ,p/)  is  the  index  set  of  all  possible  insertions  of  p  out  of  the 
r  components  of  D"  into  the  graph  jD,  divided  into  groups  with  p;  of  the  p  components 

^  factor  is  the  number  of  ways  in  which  we 

can  choose  p  components  out  of  the  r  components.  I  is  the  number  of  groups  into 
which  we  divide  these  p  components.  The  subsequent  l  —  1  combinatorial  factors  in 
(2.23)  are  the  number  of  ways  into  which  the  p  components  can  be  divided  into  l 
groups.  Components  in  each  of  these  l  groups  is  inserted  into  a  distinct  edge  of  D. 
After  insertion  the  p,-  components  can  be  permuted  in  pt!  ways,  which  accounts  for 
the  factors  pj!  •  •  -p\\.  Note  again  that  we  should  have  divided  by  the  factorials  of  the 
number  of  identical  pd  s,  but  since  we  are  interested  here  in  an  upper  bound,  we  ignore 
this  factor,  h  =  h(p,  D)  =  —c(D)d/ 2  +  a(D ,  d)  is  the  same  as  in  (2.3).  The  following 
lemma  relates  the  constant  C(|m|,  D,  e,  i )  with  the  constant  C(m,  jD,  d)  of  Lemma  2.2. 
Denote  by  I1  the  subset  of  /  that  contains  all  the  insertions  of  the  subgraphs  of  D" 
which  leave  the  number  of  twin  pairs  in  D  unchanged. 


in  the  i— th  group,  i  =  1 The 


Lemma  2.3 


I.,  ..  J  C(m,D,d )  l<d<3 

C(\m\,  \D  |,  e,  c 4  <  d  <  7  and  i  <E  If 

where  b(d ,  /),  e)  is  as  given  in  Lemma  2.2.  When  i  G  /  —  I\ 


C(|m|,|Z)|,6,i)< 


C'(\m\,D, d)  <  oo 

C{m,  D,  d)b(d,  D,  t)0{-  ln(e)) 

C(m,  D,d)b(d,D,  e)O(e) 


d  =  4,5, 
d  =  6, 
d  =  7. 


(2.24) 


(2.25) 


Proof  of  Lemma  2.3.  If  1  <  d  <  3,  then  in  this  case  there  is  no  divergence  in  e. 
So  we  can  replace  all  the  epsilon  by  zero  and  carry  out  the  calculations  as  described 
prior  to  the  statement  of  Lemma  2.3.  In  case  4  <  d  <  7,  note  that  each  twin  pair 
contributes  a  factor  b(d.  D,  e)1^^!  to  the  divergence  of  the  moment  Em{J^)h  When 
i  £  /',  the  number  of  twin  pairs  remains  unaltered  and  the  calculation  as  above  gives 
the  constant  as  in  the  statement  of  Lemma  2.3.  When  i  (E  I  —  I1,  the  number  of  twin 
pairs  in  the  graph  is  reduced.  When  a  twin  pair  is  disturbed  due  to  the  insertion  of 
some  of  the  subgraphs  of  D",  the  integrals  associated  with  this  twin  pair  will  diverge 


16 


at  a  slower  rate.  To  see  this,  consider  the  case  when  a  component  of  D"  is  inserted 
into  the  edge  leading  to  the  exit  node  labelled  U,  in  a  twin  pair  labelled  t\  and  £2- 
The  twin  pair  which  was  contributing  to  the  divergence  of  the  moment  at  the  rate 
e-(d-4)/ 2  ^see  2.9),  now  gives  a  term  of  the  form 
1  pT  pT  ptjAt 

(2^yd/2  J  J  dhk-idt2H{\ t2k-i2k-i\>t}  Js  ds-ifak-i  +  hk  ~  2s)  (2.26) 

which  is  0{—  ln(e))  for  d  =  6,  0(e-1/2)  for  d  —  7  and  stays  bounded  for  d  <  5.  It  is 
easily  seen  from  the  calculations  that  insetion  of  more  components  will  lead  to  more 
integrals  in  (2.26).  The  result  will  be  that  the  term  will  remain  bounded  in  t  for 
d  <  7.  This  proves  Lemma  2.3. 


Writing  the  result  of  the  above  Lemma  succintly,  we  can  say  that  <7(|m|,  |D|,  e,i)  < 
C(m,  D,d)b(d ,  D,  e),  for  all  i  6  I.  The  cardinality  of  /  is  |D|(|D|  —  1)  •  •  •  (]D|  —  /+  1). 
Using  these  two  facts  and  simplifying,  we  can  bound  (2.23)  by 


Dev2j  D"evk  p= 1 


V 


p!x 


n(  p,  nif  1  n  -  1)  •  •  •  (|T)|  -  /  +  1)  ~  k  . 

xEE  C(m,D,d)Hd,D,e)T  {h+ l)(h +  2)...(h +p)  Id" 


(2.27) 

1 

|m|r_p  ’ 


where  C,  b  and  h  are  as  in  Lemma  2.2. 


We  now  derive  a  bound  for  h.  Our  aim  is  to  show  that  the  fraction  |jD|(|T>|  — 
1)  •  •  •  (|jD|  —  /  +  1  )/h(h  +  1) . . .  (h  +  p),  appearing  in  (2.27)  is  bounded  by  a  constant 
depending  only  on  k. 


Claim:  j  <  h{d,D)  <  \D\. 


Proof  of  Claim:  Since  d  <  8,  h  >  —4 c(D)  +  a(D,d )  =  —4 c(D)  +  \D\  —  2 k(D). 
Obviously  h(d,D)  <  \D\.  k{D )  equals  zero  when  1  <  d  <  3,  and  equals  the  number 
of  twin  pairs  for  4  <  d  <  7.  We  divide  the  2 j  exit  nodes  in  D  into  three  categories: 
The  ones  that  form  twin  pairs,  the  ones  that  contribute  to  c(D),  and  the  rest.  For 
each  twin  pair,  we  can  associate  three  distinct  edges  in  D.  With  each  group  of  nodes 
that  contribute  one  unit  to  c(D ),  we  can  associate  at  least  six  edges  (for  it  takes  at 
least  four  exit  nodes  with  no  twin  pairs  to  add  to  c(D)).  With  the  remaining  exit 
nodes  we  can  associate  at  least  one  distinct  edge.  Thus  we  have, 

\D\  >  6c(D)  +  2k(D)  +  (2j  -  4c(D)  -  2k(D)) 

>  2j  +  2c(D)  +  k(D) 


So, 


h  >  -4 c(D)  +  \D\  -  2k(D)  >  2 j  -  2 c{D)  -  k(D) 

(j  -  HD)) 

>  2j  —  2— — — ^ — —  - 


HD)  >  i, 
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where  we  used  the  fact  that  c(D )  <  ( j  —  k(D))/2  (see  proof  of  Lemma  2.2). 

Using  the  above  claim  and  the  fact  that  \D\  <  4 j,  we  get  the  bound  |D|(|D|  — 

1)  •  •  •  (\D\  —  l  +  l)/li(h  -f  1) . . .  (h  +  p)  <  4fc,  since  l  <  p  <  r  <  fc.  This  together  with 

the  inequalities  p!  <  kl,  | Alp.  /)|  <  (  ^  j  )  in  (2.27)  yields  the  bound 

Em(l,Vl)k(jy  <  Em{l,%)kEm(jy  +(4C,)‘  H  £  X)  C(m,D,d)b(d,D,e)Tkx 

Dev2j  D"evk 

x  E(  p  )  E  (  111  )  Id”,  (2.28) 

where  C,  =  (1  V  |m|-')(l  V  T).  Since  T)*L,  (  (  j  <  2‘,  we  get 

<  Em{i,%)kEm(jy  +  (16^  k\  Em(i,p,)kEm(jy 

=  (1  +  Ck  k\)Em(l,  m)kE„(jy.  (2.29) 

^,From  (2.16),  (2.17)  and  (2.29),  we  conclude  that 

EZj(f,n,)k  <  ||/||*(1  +  k\Ck)Em( l.i,,)*.  (2.30) 

Tightness  of  the  marginals  now  follows  from  the  boundedness  of  the  moments. 

To  complete  the  proof,  we  show  that  for  any  positive  function  /  £  such  that 

ll/ll  +  ||A/||  <1, 

£££[</,*>  -  (f,v,)}t  <  CT(t  -  ‘Y-  (2.3i) 

We  know  that  (see  Dawson  (1993),  Theorem  4.7.2,  Proposition  7.3.1) 

Em[(f,pt)-U,r,,)}k<CT(t-s)\  (2.32) 

Note  that  we  cannot  apply  the  martingale  approach  of  Proposition  7.3.1  from  Dawson 
(1993),  since  the  measure  M  appearing  there  is  no  longer  a  martingale  under  the  new 
law.  We  can,  however,  follow  the  method  employed  in  proving  Theorem  4.7.2  there, 
which  involves  expanding  the  fourth  power  of  the  difference  inside  the  expectation  and 
then  carrying  out  a  direct  computation,  using  the  boundedness  of  /,  the  contraction 
property  of  the  Brownian  motion  semigroup  and  the  relation 

II TJ  -  T,f\\  ss  ||A/|1((  -  s).  (2.33) 

The  calculation  essentially  involves  a  clever  juggling  of  the  ranges  of  integration  of 
the  time  variable  associated  with  the  exit  nodes.  To  analyse  (2.31),  expand  the  fourth 
power  of  the  difference  on  the  left  after  writing  the  expectation  E  as  a  ratio.  Each 
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term  in  the  numerator  involves  evaluating  a  moment  of  order  2j  -f  4.  We  divide  the 
graphs  contributing  to  the  2 j  -t-  4  moments  into  two  parts  as  done  while  evaluating 
Em{  1, The  Par^  consisfs  °f  attaching  separate  graphs  from  V4  to  graphs 
from  V2 j.  This  will  give  a  term  CT(t  -  s)2.  The  second  part  consists  of  those  graphs 
that  can  be  constructed  by  inserting  graphs  from  V4  into  the  edges  of  graphs  from 
V2j  as  done  earlier  in  this  proof.  Consider  the  terms  in  the  expansion  that  correspond 
to  the  same  graph  from  V2 j  into  which  a  graph  from  V4  (this  graph  will  be  same  for 
all  the  terms  except  that  the  time  variable  corresponding  to  the  exit  nodes  will  be 
different)  is  inserted.  The  factor  coming  from  the  graph  from  V2j  is  common  to  all 
the  terms  and  can  be  factored  out.  The  inserted  parts  can  now  be  analysed  exactly 
as  in  the  proof  of  Theorem  4.7.2  (Dawson,  1993).  This  gives  us  a  factor  ( t  —  s)2.  The 
remaining  terms  can  evaluated  exactly  as  in  the  proof  of  Theorem  1.2.  This  will  give 
us  a  bound  for  the  required  expectation  of  the  form 

{CT  +  C(T,  d,  D,  |m|))(t  -  s)2  :=  C'T{t  -  s)2, 

which  completes  the  proof  of  tightness.  ® 


Proof  of  Theorem  1.3.  a.  So  as  to  avoid  some  cumbersome  integrals  involving 
logarithms  in  even  dimensions,  we  shall  prove  the  result  only  for  odd  values  of  d  <  7. 
Filling  in  the  even  values  requires  patience,  but  no  further  ingenuity.  We  start  with 

/T  rti 

dti  J ^  dt2(6(x  -  y),r)tl(dx)T}t2(dy)}, 

and  first  evaluate  Em(Jc). 

Em(Jc)  =  2j^  dti  J  dt2  J  m(dx1)m(dx2)ptl(x1,z)pt2(x2iz)dz+ 

2  J  m(dx)ps(x,y-l)ptl-s(yi,  z)pt2-s{yi,  z)  (2.34) 

The  first  term  can  be  written  as 

2Jc  dtx  L  (2-35) 

where  we  denote  by  /  the  density  of  m  w.r.t.  the  Lebesgue  measure.  The  second 
term  gives 


tT  rt\  ft2 

4  /  dt\  /  dt2  /  ds\m\(ti  +  t2  —  2s) 
Jt  J  0  Jo 


- d/2 


which  equals 

8|m| 

(<J-2)(d-47 


^  - (2r  -  + (T ~ e)e_ 


d—A 

2 


.  (2.36) 
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hspace*1.0truein  (v)  (vi)  (vii) 

Figure  4:  A  graph  for  the  compuation  of  Em(l,r]t)J £. 


We  now  proceed  to  evaluate  Diagramatically,  we  can  express  this  moment 

as  in  Figure  4.  The  first  two  graphs  add  to  give  \m\E(Jc).  The  third  graph  corresponds 
to 

rT  rt\-L  r 

4  j  dti  J ^  dt2  J 1  ds  Jm(dx)m(dy)ps(x,yi)pL-s(yl,w)ptl_s(y1,z)pt2(y,z)dyidzdw 

tT  rti-t  dti  At)  r 

=  4  /  dti  dt2  /  ds  m(dx\)m(dx2)pt-\+t2{x\,x2)  :=  Ae2(t,T)  (2.37) 
Jt  Jo  Jo  J 

Similarly,  the  fourth  graph  gives  a  contribution  of  A\(t,T).  The  fifth  graph  gives  a 
contribution  of 


(d  —  2)(d  —  4) 


_fci  r(T-e)2 
£  !  [^“ 


(T-t-tf 


(d  ~  6) 


r  d—6 

\-3e~— + 


{l(2T  -  e)-"  -  (2 T  -2 1-  e)~^]  +  8[(T  -  -  r'^]}}  J  (2.38) 

The  remaining  terms  can  be  similarly  evaluated.  The  point  to  note  here  is  that  for 
4  <  d  <  7,  the  only  terms  that  will  matter  in  the  e  — *  0  limit  are  those  of  order 
e-(c2-4)/2  terms  we  have  not  evaluated  explicitly  here  will  diverge  at  a 

slower  rate.  In  the  case  d  <  3,  all  terms  containing  multiples  of  t  will  vanish  leaving 
a  ratio  that  will  look  like  the  one  given  in  the  statement  of  the  theorem.  In  this 
case  terms  coming  from  all  the  graphs  will  matter.  Furthermore,  since  each  of  the 
integrals  we  evaluated  is  positive  for  all  e,  the  result  follows. 

b.  To  prove  part  (b),  we  will  write  a  finer  version  of  (2.23)  for  the  case  k  =  1.  See 
proof  of  part  (a)  above  for  illustrations.  Fix  n  >  1.  As  e  — ►  0,  only  the  terms  with  the 
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highest  order 
we  can  write 


of  divergence  in  e  will  be  relevant.  From  Lemma.  2.2  and 

Em[(l,m)(JT] 


(2.23),  (2.27), 
(2.39) 


as  e  — >  0.  As  is  now  routine,  we  divide  the  numerator  into  a  sum  where  a  graph  in 
V 2n+i  is  obtained  by  considering  the  graphs  in  V 2n  and  adding  to  them  the  graph 
in  T>i  as  a  seperate  component.  This  on  dividing  by  the  denominator  gives  the  term 
=  \m\-  We  now  focus  on  the  second  part.  In  course  of  the  derivation  of 
(2.23),  if  we  use  the  fact  that  the  effect  of  insertion  also  restricts  the  range  of  the 
time  variables  in  the  graph  D,  (see  argument  following  (2.22)),  then  we  would  have 
ended  with  a  factor  Th+P  —  (T  —  t)h+p  instead  of  Th+P.  Further  as  e  — *•  0,  only  those 
graphs  in  V2n ,  matter,  for  which  there  are  n  twin  pairs  (see  Lemmas  2.2, 2.3).  We 
index  the  set  of  such  graphs  by  £2n-  For  any  D  6  £2ri,  k(D)  =  n,  and  C(D)  =  0.  So 
h(D )  =  |D|  —  2 n.  Now  for  the  combinatorics  part.  The  number  of  ways  in  which  we 
can  insert  the  only  edge  of  the  graph  D"  (E  V\  into  a  graph  D  €  £2n  so  as  not  to 
disturb  the  number  of  twin  pairs  (see  Lemma  2.3)  is  |jD|  —  2 n.  So,  in  the  limit  e  — >  0, 
(2.39)  behaves  like 


\m\  + 


m 


-•  C(m,D, <QK4  D, e)(T™  -  (T  -  t)h*') 

T )ncf„_  C(m,D,d)b(d,D,e)Th 


(2.40) 


As  n  -*  oo,  h  =  (|D|  -  2n)  — »  oo,  since  3n  <  \D\  <  4 n,  and  so  (2.40)  converges  to 
\m\  +  |m|_1r. 

This  proves  Theorem  1.3.  ® 


Proof  of  Theorem 

we  can  write 


1.4.  Using  the  same  arguments  as  in  the  proof  of  Theorem  1.3, 


*t<  UmY 


^[(i,i.)V<)“)  t>1 

EJJT  ’  ~ 


(2.41) 


as  e  — *  0.  Again  proceeding  as  in  Theorem  1.3,  we  rewrite  (2.23),  by  using  the  precise 
limits  of  integration  while  integrating  the  time  variables.  Then  we  simplify  to  get 
a  form  similar  to  (2.27.  As  above  we  note  that  as  e  — >  0,  only  graphs  in  £2n  will 
matter.  Now,  in  the  n  — »  oo  limit,  since  h  — *  oo,  the  only  combinatorial  factor  that 
will  balance  with  the  denomintor  (see  (2.27)  are  the  ones  for  which  l  =  p,  that  is 
Pi  =  1,  for  i  —  1, . . .  ,p,  i.e. 


m\p\- i)---(i/?i-*+i)  f  i  i fi=p  f242) 

(|T»|  +  l)(|D|+2)---(|D|+p)  1  0  otherwise 

So,  in  the  limiting  procedure  prescribed  in  the  statement  of  the  theorem,  (2.41) 
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behaves  like 


Em{  1)  Vt)k d 


l^D"evk  Bm(l,Vt)k\D"  LJ.=i  />(r,i,p,D")j^j~SDe£2nC(m,  £>,  d)b(d,  D.  e)Th 


ZD&*.C(m,D,J)b(i,D,')T* 


(2.43) 


where  P(T,t,p,D")  is  a  polynomial  of  degree  p  in  t,  T.  Thus  we  have  proved  that 

eJ-U  7-4  // , — n  _ t  ?  n  1 


D"evk 


v- 1 


(2.44) 

Repeating  the  above  arguments,  we  can  show  that  the  mixed  moments  of  the  process 
Xt  are  unique  under  the  given  limits.  This  proves  Theorem  1.4.  ■ 


Proof  of  Theorem  1.5.  It  is  clear  that  exponential  integrability  will  follow  once  we 
get  an  appropriate  bound  on  the  the  moments  of  (/,  Tjt)  k  <  ||/||oo(l,f?t)fc-  From  (2.23), 
(2.27)  and  the  arguments  in  Theorem  1.4,  we  have  P(T,t,p,  D")  <  P(T,T,p,D")  = 
Tv  <  (1  VT)fc.  P(T,t,p,D")j^-  <  c*,  where  Cx  =  (1  VT)(1  V  \rn\~1).  This  together 
with  (2.44)  gives 

lim  \imE^r(l,rit)k  <  Em{l,r]t)k{\  +  (2 Cr)k)  <  CkEm{l,Tjt)k .  (2.45) 

n—+oo 

This  implies  that 

Jim  lim EZ#{f,%)k  <  (C||/||)‘£„<1,„)‘  (2.46) 

The  result  now  follows  from  the  exponential  integrability  of  the  Feller  diffusion;  i.e. 
from  the  fact  that  Eme0<'1'Vt'>  exists  and  is  finite  for  all  t  <  0_1. 

Proof  of  Theorem  1.6.  By  definition, 


EXT 

x,T,a 


)  = 


£x,r,q(e-c*t+A^) 

Ex,T,a(eX,lT) 


Let  v 


a, A  ,c 
t 


denote  the  unique  solution  of  the  following  differential  equation: 


dvt  2 

—  =  avt  -  v  -  A  ;  v0  = 
at 


c. 


Iscoe  (1986)  has  shown  that 

Em,T.Ae~{f  ’Vt)+X  ti{4’'rls)ds)  = 


(2.47) 


(2.48) 


(2.49) 
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where  is  the  unique  solution  of  the  non-linear  initial  value  problem. 

—  =  (A  +  a)ut  -  u2  -  \<j> 

St 

Uo  =  /• 


(2.50) 


When  the  functions  /  and  <f>  are  constants,  then  u  no  longer  has  a  spatial  depen¬ 
dence.  The  denominator  in  (2.47)  can  be  written  as  exp(-sv£’  ’  ).  Using  the  Ma.rkov 
property  and  homogeneity  the  numerator  can  be  seen  to  be 


ex-p(-xvt 


\  X  <2,X,0 
a1Aic-tvT_t 


)• 


(2.51) 


Proving  Theorem  1.6  now  involves  verifying  that  the  denominator  and  the  numera¬ 
tor  in  the  various  cases  satisfy  the  requisite  differential  equation.  For  example,  the 
denominator  in  case  (»)  on  the  right  can  be  easily  seen  to  be  exp(-xuT'  ’  )•  We  leave 
the  straightforward  verification  of  the  other  cases  to  the  enthusiastic  reader  and  s^ 
conclude  the  proof  of  Theorem  1.6. 


Proof  of  Corollary  1.7.  The  corollary  follows  immediately  by  noting  that  m  each 
case,  as  T  approaches  the  critical  value,  the  denomiator  tends  to  +oo,  while  the  nm 
merator  stays  bounded. 


3  Some  results  for  the  reinforced  Feller  diffusion 


Any  number  of  moment  formulae  and  relationships  can  be  gleaned  from  the  formulae 
in  Theorem  1.6.  Here  are  a  few  illustrative  examples. 


Corollary  3.1  Let  7  be  as  in  Theorem  1.6.  Let  A  >  0. 

(1) 


(vi) 

Var^Tfi(Xt) 
(lii)  If  7  =  0,  then 


x  _  g  sec2(\/AT) 

Ex,Tfi\Xt)  -  sec2  (^(T-t))' 

2x  se^fV A  f)[tan(y/A  T)  -  tan(\A(T  -  t))] 
=  71  sec4(\/A(T  —  t)) 

(1  _  2 (T-t))x 

r?  A  MM  _  1 _ - — — 

•e,x,T,cA'AU  —  ^  _  a 


(3,1) 


(3.2) 


(3.3) 
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(iv)  7/7  =  0,  and  a  <  0,  then 

Jim  Ex\rJe-'*‘)  =  £.(«-*•)  =  (c 

where  the  right  side  of  the  equality  is  the  Laplace  transform  of  a  Feller  diffusion. 

(v)  Let  a  =  —2y/\.  Then. 

lim  E~Tfi(Xt)  =  Ex.ooAXt)  =  xeat.  (J 

T— ►  00  5  * 


Remark:  Result  (w)  says  that,  depending  on  A  and  a,  the  expected  reduction  in 
mass  in  time  for  the  reinforced  Feller  diffusion  and  a  simple  Feller  diffusion  with  a 
mass  annihilation  term  are  almost  identical  for  large  T.  The  distributions  in  the  two 
cases  are,  however,  quite  different. 

Proof  of  Corollary  3.1.  The  corollary  is  an  immediate  consequence  of  Theorem  1.6 
and  the  equalities 

=  -fjLJ*-**')  1=0  (3.6) 


£*,7>a'<2  =  U° 


a  f^-cXt 


Corollary  3.2  (Extinction  Probabilities)  Let  A  >  0  and  let  7  be  as  above. 

0) 


Px,T,a(Xt  =  0)  =  exp 


1  —  e° 


(ii)  If  y  >  0  and  0  <  T  <  ^  (f  —  tan  *(2^7))  >  then 

Px,T,a(Xt  =  °)  =  exP{-\/Mtan(7f  -  V^7)  +  tan(tan_1(-^)  +  y/yT)]}.  (3.9) 


(in)  If  7  =  0  and  0  <  T  <  2/a  if  a  >  0,  then 


p*,TJXt  =  0)  = 


(iv)  If  7  <  0  let  (3  =  y/~T-  Then 


Px,T,a(Xt=°)  =  exPi-^ 


1  -  c«‘  ifi  -  |)e«T  -  (/?  + 1)- 
1  -  e»‘  (!)  -  f)e®T  +  (0  +  f ) 


(3.10) 


(3.11) 
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Proof  of  Corollary  3.2  The  corollary  follows  from  Theorem  1.6  and  the  equality 


(3.12) 


4  Appendix 

For  completeness,  we  shall  now  discuss  the  finiteness  of  the  moments  of  the  approxi¬ 
mate  self  intersection  local  times  J£.  While  this  result  is  similar  to  others  that  have 
appeared  in  the  literature,  it  is,  in  fact  new. 

Theorem  4.1  E(Jcy  <  oo  for  all  j  >  1,  T  <  oo  and  d  <  7. 

Proof.  We  only  give  a  sketch  of  the  proof.  The  reader  is  referred  to  Dynkin  (1988), 
Rosen  (1992)  and  Feldman  and  Iyer  (1996),  where  similar  results  can  be  found,  for 
the  details  on  evaluating  moments  similar  to  the  one  below. 

Let  -ftf(si,s2,t/i  -  2/2)  =  Is  dfi<ft2pta+*2-si-s2(2/i  -  Vi),  where  B  =  [0,T]2  n  {(fi,t2)  : 
\ti  —  t2 1  >  e}.  Recall  that 

J£  =  /  {8(x  -  y),T]tl(dx)r]t2(dy))dtldt2, 

Jb 

so  we  can  write 

(Jey  =  n  /  {S(z2i-i  -  Z2i),Vt2i-i(dz2i-i)rh2i{dz2i))dt2i-idt2i. 

»= i  Jb 

E(jy  =  @d-  Consider  a  term  Cd-  We  separate  our  analysis  of  Cd  into  the 

following  cases: 

(1)  If  the  exit  nodes  labelled  (f2,_i,  z2;-i)  and  (t2i,  z2i)  are  paired,  that  is  they  have 
an  immediate  common  ancestor  labelled  (sk,yk)i  then  integrating  w.r.t.  z2j-  and  ^2t-_  1 
and  using  the  fact  that 

I  S{z\  -  z2)f(z1)g(z2)dz1dz2  =  /  f{w)g{w)dw,  (4.1) 

J$2d  Jfcd 

we  obtain  a  factor 

i?f(-Sfc,-s/c,0)  =  I  dt2i^dt2ipt2i_1+i2i-2sk  <  c\B\e~d/2.  (4.2) 

J  B 
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So  we  can  bound  the  contribution  coming  from  all  the  paired  nodes. 

(2)  If  at  least  one  of  the  exit  nodes  £21-1  or  z3i  is  connected  to  a.  root  (entrance)  node 
labelled  (0, £/),  then  we  get  a  contribution  /  m(dx )  JB  dt2i-idt2iPt2i_1+t2i-Sk(yk  —  %i)- 
Using  the  bounded  density  of  m  we  integrate  out  xi  to  see  that  the  above  term  is 
finite. 

(3)  Now,  of  the  remaining  exit  nodes  whose  space  variables  z? s  a.re  left  to  be  inte¬ 
grated,  none  of  them  are  paired  or  are  connected  to  a  root  node.  (These  have  been 
taken  care  of  in  the  above  two  cases).  Suppose  there  are  2 r  remaining  2,-’s.  If  z2i  and 
Z2%-\  are  connected  to  the  internal  nodes  labelled  (sk,yk)  and  (si,yi),  then  integrating 
with  respect  to  z2i  and  z2i- 1  using  (4.1)  gives  a  factor  R(sk,  si, yk  —  yi).  Since  there 
are  2 r  zds  to  be  integrated,  we  will  get  r  factors  that  look  like  R(sk,si, yk  —  yi).  Each 
(s k,  yk)  occurs  in  exactly  2  of  the  R's.  Make  suitable  change  of  variables  (see  Feldman 
and  Iyer,  1996,  for  details)  w{  =  yk  -  yi,  such  that  {y{l , . . . ,  yip,  , . . . ,  w2r-P}  spans 
di2rd.  Since  each  function  R  in  the  integrand  has  an  argument  that  is  a  difference  of 
the  yd s,  upon  change  of  variables  all  the  R’s  will  have  arguments  depending  only  on 
the  w’s.  Furthermore,  each  w  occurs  in  exactly  two  of  the  R’s.  Also  to  be  noted  is 
the  fact  that  one  cannot  have  W{  —  wj  occuring  as  the  sole  argument  in  two  of  the 
R’s. 

First  integrate  out  all  the  variables  , . . . ,  y*p.  Now  using  the  bounded  density  of  m, 
if  necessary,  one  can  integrate  out  all  the  internal  space  variables.  This  will  leave  us 
with  the  factors  R  and  the  time  integrals  only.  Now  we  integrate  with  respect  to  the 
w’s.  If  one  has  a  factor  /  R2(si,  S2,  W\)dwidsids2,  it  can  be  written  as 

[  ds\ds2  f  P«i+t2+t3+t4-2si-2s2(0)  H  dti  <  T2  f  (ti  +  ^2  +  ^3  +  ^4)  d/2  IJ  dti 
J  Jb 2  i= 1  ^toT)4  ~i 

<  00  for  d  <  7.  (4.3) 

A  factor  of  the  form 

J  dsids2dsz  J  R(si,S2,Wi)R(s-i,s3,W2)R(s2,s3,W2  —  w-i)dwidw2dw3  (4.4) 

arising  out  of  the  configuration  shown  in  Figure  5  is  bounded  by  T3/j0Tj6(ti  +  . . .  + 
te)~d/2dti  . .  ,dte  is  finite  for  d  <  11.  Similarly  we  can  check  that  the  more  one  has 
to  integrate  w.r.t.  the  w’s,  the  smoother  the  integral  becomes,  the  worst  case  being 
given  by  (4.3).  We  conclude  the  proof  of  Theorem  4.1.  H 
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Figure  5:  Computing  the  factor  (4.4). 
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